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Abstract 

We investigate the Green's functions G(x;x') of some second order 
differential operators on R d+1 with singular coefficients depending only 
on one coordinate xo- We express the Green's functions by means of the 
Brownian motion. Applying probabilistic methods we prove that when 
x = (0, x) and x' — (0, x') (here xq = 0) lie on the singular hyperplanes 
then G(0, x; 0, x') is more regular than the Green's function of operators 
with regular coefficients. 



1 Introduction 

We discuss Green's functions of some second order differential operators with 
singular coefficients appearing in quantum physics. 

As a first example consider the Lagrangian for a scalar field in (d + 1)- 
dimensions interacting with gravity 

C = g^ v d^d v (t> + (m 2 + £R)cj) 2 (1) 

where is the metric tensor and R is the scalar curvature. Such a Lagrangian 
with m = and the minimal coupling £ = appears also in the theory of 
structure formation (cosmological perturbations) [1]. We discuss the Euclidean 
version of a spatially homogeneous metric (we write x = (t,x) or x = (xo,x) 
depending on whether the first coordinate has an interpretation of time or space) 

ds 2 = dt 2 + g ]k (t)dx j dx k 

The Laplace-Beltrami operator resulting from the bilinear form in eq.(l) reads 

A s = \g-*dr<jrgtd v ) (2) 
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(here g = det(g fil/ )). In cosmological models gjk — t 2a and gi k ~ t~ 2a when 
t — > with a > 0. Such a singular behavior can appear also in models describing 
collapse phenomena in general relativity [2]. 

As a second example we consider quantum mechanics on a (topologically 
trivial) manifold with the Hamiltonian 

H = -A g + U(x ) (3) 

(in some global coordinates x = (xo,x)). 

The Green functions of (Euclidean) quantum scalar fields (1) with m = 
and the minimal coupling £ = are solutions of the equation 

-A g G = g-h (4) 

These Green functions are also relevant for classical field theory because they 
describe a propagation of disturbances. In quantum mechanics (3) we are in- 
terested in the propagator kernels 

exp(-TH)(x ,x;x' ,x') (5) 

where t is purely imaginary. 

In this paper we prove that if the coefficients of the Laplacc-Bcltrami op- 
erator have a power-law singularity at a certain point t = to then the Green 
functions G(to, x; to, x') are more regular than the ones of operators with reg- 
ular coefficients (for regular coefficients the Green function can be expressed 
by the geodesic distance [3] [4]). In quantum field theory these Green functions 
have the meaning of expectation values of quantum fields at equal times. In 
quantum mechanics the propagator (5) will have an anomalous behavior in r. 
The Green function (4) can be obtained from the propagator (5) by means of 
an integration over r. 

2 The Green's functions 

Let us change coordinates 

The Laplace-Beltrami operator (2) takes the form 

2A g = g- 1 d 2 ri +g^d j d k (7) 

The bilinear form in eq.(l) determines an operator A which is of the same form 
as H in quantum mechanics (eq.(3)) 

A=-A g + w (8) 
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Hero, w = jm 2 + for the scalar field and w = U for quantum mechanics. 
The Green's function of A is a solution of the equation 

-{dl + gg ]k d 3 d k - W)G = 2S(r] - »/)<f(x - x') (9) 

where we write 

W = gw (10) 
Together with eq.(9) we consider the differential equation 

-d T P T =AP T (11) 

with the initial condition Po(i], x; ?/, x') = 5(rj — rj')S(x — x'). Eq.(ll) defines 
the transition function of a stochastic process [5]. 

We can formulate the problem of solving the equation 

AG = 5 (12) 

as a problem in the Hilbert space of square integrable functions L 2 (dr]dx) [6]. We 
assume that W is a non-negative function. The operator A can be considered as 
a self-adjoint non-negative operator in L 2 if ggi k and W are locally integrable 
functions (then we can define the Friedrichs extension [6] of the symmetric 
differential operator (8)). The transition function P T of eq.(ll) can be defined 
as the integral kernel of exp(— tA). Then, the kernel of the inverse 



A- 1 



f-OO 

I drcxp(— tA) 
Jo 



is the solution of cq.(9). It follows that the Fourier transform G of G has the 
representation 

/•OO 

G( V , V ';p)= / dTP T (v,r/,p) (13) 
Jo 

where P is a solution of the equation 

-d T P T = AP T (14) 

with the initial condition Po(r],r]',p) = 5(r] — 77') ( the fundamental solution). 
Here 

A = -\d 2 n + \ Pj gi k {Ti)g{ri)Pk + W = -\d 2 r) + V{rf) + W{rj) (15) 

Eq.(14) is a Schrodinger-type equation with the Hamiltonian A and the potential 
V + W where 

v(v) = lpjg Jh (v)g(v)pk = pVp(v) (16) 

If the potentials V and W belong to Lj oc (dr]) then A is a well-defined essentially 
self-adjoint operator in L 2 (drj)[7]. 
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We can express the kernel P by means of the Brownian motion b (the 
Feynman-Kac formula [8]; a discussion of the probabilistic representation for 
singular potentials can be found in [9]) 



p t (v,v',p)= nri 

E[S (rf- V -b (r)) exp (- f Q V ( V + b (*)) ds - f Q W (r, + b («)) da)] 1 '> 

where E[.] denotes an average over the Brownian paths. Now, the kernel of 
exp(— tA) has the representation 



P T (77,x,r/,x') = (2TT)- d J dpexp(ip(x'-x)) 

E[S (v'-V-b (r)) exp (- f Q r V fa + & (*)) ds -J^W(r] + b (a)) ds)} 



(18) 



In order to eliminate the 5 function in eq.(17) it is useful to express the expec- 
tation value over the Brownian motion by means of an expectation value over 
the Brownian bridge 7. Let q be a path connecting 77 with rf 

q(-) = v + (v' -v) 1 + V^i(-) (19) 

T T T 

where 7 is the Gaussian process on the interval [0, 1] (the Brownian bridge) 
starting from and ending in with the covariance 

Eh(*h(8')] = Ai - s) 

for s' < s. Then, eq.(17) can be rewritten in the form [8] 

Pr(v,v',P) = (27rr)-5 C xp(-^:(?7' - r]) 2 ) 



£[exp (-T £ ds (V (q («)) + W(q (s)))) 



(20) 



Applying the Jensen inequality ( see [10]- [11]) to the E[..] integral we obtain 
the inequality 

> (2nT)-i exp(-£(i/ - V ) 2 ) (21) 
exp (-r £ dsE[V (q (a)) + W(q («))]) = P L 

This integral is 
P?(V,r/,P) = exp(-£(i/ - V ) 2 ) 

exp (-r d a /dj/(2™(l - ex p(- 5 ^ y )(F + (r? + s(r/ - r?) + 

(22) 

As a simple application of the inequality (21) we note that if 

V + W < A'p 2 + B' (23) 
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then 



PHv,rf,p) > (2^r)-icxp(-JL(V-ry) 2 )cxp(-rA'p 2 -r J B') (24) 

Hence, we obtain a bound from below by the transition function for the d- 
dimensional Brownian motion. 

On the other hand we may apply the Jensen inequality in the opposite 
direction to the s-integral 



< (2irr)-i cxp(-^(V - ry) 2 ) £ dsE[c X p (-rV (q (s)) - tW (q (a)))] = P u 

(25) 

This integral takes the form 
= (2^T)-^exp(-27(77'-??) 2 ) 

Ji ds J dy(27rs(l - s))~ L 2 exp(- 5 ^) cxp (-r(V + W) (r, + a(rf - rj) + y/ty)) 

(26) 

If 

V + W > Ap 2 + B 

then 

< (2^r)-^exp(-^(r/'-?7) 2 )cxp(-TAp 2 -TS) 1 ' 

Hence, we estimate the transition function from above by the Wiener transition 
function. 



3 Scale invariant metrics 

We consider in this section a power-law cosmological expansion. Such an ex- 
pansion is an exact solution of coupled Einstein equations for a metric and for 
the scalar field with an exponential self-interaction . Some consequences for a 
structure formation with such an expansion are discussed in [12] [13]. If gjk(t) 
has an isotropic power-law behavior then V is scale invariant. Let us assume 
here that V and W are nonnegative and scale invariant around rj — (there is 
nothing special in the choice of rj = as a singular point, see a discussion at 
eq.(44)) 

V jk (Xrf) = \ 2 "V ik (ri) (28) 

and 

W(\ri) = X 2a W(r]) (29) 
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Let us denote 9 = r 2 r\. We apply the scaling properties of the Brownian 
bridge (19). Then, for V of the form (28) and W (29) we obtain 

Mv, V',P) = (27rr)-5 eM-^W - V?) 

£[exp (-r^ Jl pVp (6 + s(9> - 6) + 7 (*)) da) ( 30 ) 
cxp (-r 1+ff /J (6» + s(0' - 0) + 7 («)) ds)] 

The bounds (22) and (26) become simple if 77 = ry' = 0. Then, the bound (22) 
reads 

P T i (0,0,p) = (2 7 rr)-i 

exp ( - r 1 *" Jo / dj/(2™(l - s))-W{y) exp^^f^) 

-t^° & ds S dy(2ns(l - a))-^W(y) exp(- 5 ^ y )) 

= (27rr)-5 cxp (-T 1+ »php Si ds(s(l - s)Y - Bt 1+<7 Si ds(s(l - s)) CT ) 

(31) 

where the bilinear form h in p is defined by 

php = (2tt)-' /" rfycxp(-^)pl/(y)p (32) 

and the constant B in eq.(31) is 

B = (2tt)-* / dycxp(-^)^(y) (33) 

The integral (32) is finite if v > — | and (33) is finite if a > — |. In such a case 
the lower bound (31) is non-trivial. The upper bound (26) takes the form 

P T c/ (0,0,p)-(2^r)^/ 1 d S /dy(2^(l- S ))-i 
cx P(-^)) cx P(-^ 1+ ^(y)-^ 1+ ^(y)) 

We are interested in the Green functions (9) of the operator A which ac- 
cording to eqs.(13) and (18) are expressed by an r integration upon P T . As the 
simplest example of the integral (13) let V + W = Ap 2 + B then performing 
the r integration upon the rhs of eq.(13) we obtain 

G o (0,0;p) = (2Ap 2 + 2B)-5 (35) 

This is the standard behavior of equal-time Green's functions for the quantum 
free field. 

In eq.(31) let us first discuss the case W — B = 0. Then, the integral over 
r of eq.(31) gives the lower bound on the Green's function 

G(0,0;p)>K 1 (php)-^ (36) 
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where 



2(1 



In order to estimate the upper bound (34) (for W 
bound l^lo on V, i.e., for p 7^ 

pVp > p 2 \V\ > 



(37) 

0) let us assume a lower 



(38) 



Now, we change variables in eqs.(13) and (34) (r, y) — ► (p,u) where 

p = r|p|^|t>( 2/ )| T ^ 

u = y(s(l - s))~i 
Then, the upper bound takes the form 

G(0,0;p)<|p|- a -/ oo dp(27rp)-i 

^d S Jdu(2n)-Hs(l- S ))-^\V(u)\^eM-i)^p(-p 1+I ') 



(39) 



We can see that the integral on the rhs of eq.(39) is finite if — 1 < v < 00. 
We can summarize our results as 
Theorem 1 

Assume that W = and the potential V in eq.(15) is nonnegative and scale 
invariant with v > — | (cq.(28)). Then, the operator A is essentially self-adjoint 
and the integral kernel of exp(— tA) has the probabilistic representation (20). 
The Green's function G of eq.(9) can be defined as an integral kernel of A^ 1 . 
Assume that the potential V satisfies the lower bound (38) then the Fourier 
transform G(r],r]';p) of G(t),t)' ,x — x') at 77 = 77' = for any p satisfies the 
inequalities 

Kiiphp)-" < (5(0, 0; p) < K 2 \p\- 2 " (40) 

where h is defined in eq.(32), K\ and K 2 are some positive constants. 

For v < the Fourier transform G is decaying to zero faster than the Green 
function for operators with constant coefficients. As a consequence G is less 
singular than the one for operators with constant coefficients (see eq.(43) below). 

In the configuration space if W — then we can extract the r dependence 
from V using its scale invariance. Then, changing the integration variable in 
eq.(18) p = T _ 2( 1 + iy )k we can conclude that P has the form 



P T (r], x, rf, x') = T -^ 1+ ^ d -iF(T-^r], r~ V, t-^ 1+v \x - x')) 



(41) 



with a certain function F. Integration over r with a rescaled r = r|x — x'| 
brings the Green's function at equal time to the form 



G(?7,x,r/,x') = Ix-x'r^+^/dx-x'l^ry^x-xOlx-x'l- 1 ) (42) 
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It follows 

G(0,x,0,x') = Ix-x'r^+TTTT/ax-xOlx-x'r 1 ) (43) 

We obtain such a behavior in |x — x'| if we apply the inverse Fourier transform 
to the functions on both sides of the inequalities (40) . 

Let us note that if V is singular at rjo ^ (e.g., V ~ \rf — f?o| 2ly ) then all our 
results concerning the transition functions and Green functions still hold true 
but instead of setting rj = rf = we set rj = r)' = rjo (this conclusion follows 
directly from eq.(30)). So, e.g., the formula (43) reads 

G( Vo , x, r? ,x') = |x - x'|- d+ ^/((x - x')|x - x'p 1 ) (44) 

We admit now W ^ 
Theorem 2 

Let W > be scale invariant (eq.(29)) and a > — \ then (under the assump- 
tions of Theorem 1 concerning V) for any A > if |p| > A then there exist 
positive constants K\ and K 2 such that the inequalities (40) hold true. 

Proof: setting W = in eq.(34) we obtain the upper bound (39). For the 
lower bound we note that the exponential in eq.(31) is dominated by the term 
quadratic in the momenta. We change the integration variable in eqs.(13) and 
(34) 

t = r(php) i+" 

Then, we can see that for any A > if |p| > A then there exists a constant C 
such that in the exponential of eq.(31) Bt 1+ " < Cr 1+rT . Then 

f drP? (0,0,p)> 

(ph 2 p)^ J dr(2irr)-i cxp (-r 1+,y ds(s(l - s)) v - Cr 1+a ds(s(l - s)) ff ) 

From this lower bound and from the upper bound (39) we obtain the results of 
the theorem. 

If W > then the lower bound in cq.(40) cannot be true for arbitrarily small 
p because as follows from eq.(34) (V = for p = 0) 

G(0, 0, 0) < f °° dT(2n T )-i J, 1 ds J dy(2™(l - a))"* 
exp(-2^ T )exp(-T 1 +^(y)) <^ 1 ' 

If we imposed the condition that t > (which is quite artificial in the Eu- 
clidean framework) then we would need to impose boundary conditions at r/ = 
on the Brownian motion in the path integral (17). The Dirichlet boundary con- 
ditions can easily be imposed in the functional integration framework. We just 
insert the characteristic function of the positive real axis in the path integral 
(17) rejecting all the Brownian paths which leave the positive real axis. With 
the Dirichlet boundary conditions our estimates on the upper bound remain 
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unchanged whereas the estimates on the lower bound require some minor mod- 
ifications. 

Let us consider an example of a threedimensional space. By a change of 
coordinates we can diagonalize the matrix (gjk) 



9jk = SjkOj (46) 



Let a = (010203)3 and 



Sj = a,j a 2 d ri a,j 
S = a~ 3 d v a 

3<k 

Then, in the potential W of eq.(9)[14] 

gR = 6a 4 (a- 2 d ri S + 5 2 + Q) (47) 

and 



m 2 g = m 2 a 6 (48) 

id W if aj are seal 
the simplest case when all a 3 are equal, t <E R and 



We obtain a scale invariant V and W if aj are scale invariant. Let us consider 



a(t) = \t\ a (49) 

We have 

V = (I - 3a)-H\t\- 3a 

Note that for a > | the point t = corresponds to r] = —00 and t = 00 to 
?7 = 0. 
Then 

V(y) = K p 2 \y\ 2 » 
where k > is a certain constant and 

v = 2a(l -3a)" 1 (50) 

For a scale invariant metric 

W = m 2 g(r,) + £gR = dm 2 ^^ + £C 2?r 2 

De Sitter space can be obtained as a limit a — > 00. Then, we have F(r?) = 
cp 2 \r]\~3 and m 2 g = c'r]~ 2 , hence VT(^) = cr}~ 2 . This is a singular perturba- 
tion which goes beyond our analysis. It can be treated by means of the path 
integral methods. However, in such a case W needs a regularization, then a 
renormalization and a subsequent removal of the regularization [9]. The ?]~ 2 
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singularity comes also from the term gR. Hence, the results of this section ap- 
ply only to £ = 0. Then, in cq.(29) a = 3a(l - 3a)" 1 . B in eq.(33) is finite if 



\a\ < §. 

In quantum mechanics xq is interpreted as a space variable. The metric 
takes the form (d + 1 = 3) 

ds 2 = dxl + \x \ 2a (dxl+dxl) 

Then, r\ = (1 — 2a)~ 1 x \x \~ 2a ■ The Hamiltonian (3) is symmetric in L 2 {^fgdx). 
The change of coordinates x — > r\ associates with H the operator A = gH which 
is symmetric in L 2 (di]dx) 

A = —d 2 + V + W 

where 

V( V ) = C 1 p 2 \r 1 \^ (51) 

with p 2 = p\ + p\ and 

W = gU{ri) = CM^U{ri) (52) 

The anomalous behavior of P T has as a consequence 
Corollary 3 

Let P T (r),r]',p) be the fundamental solution of eq.(14) with W — and V 
defined in eq.(51). If v = l ^ 2a > — 5 then for any r > 



y rfxP r (0,x,0,x')|x-x'| 2 = (-A p )P r (0,0,p)| p=o = B\T^ 



(53) 



and 



J 'dxd»/i > r(0,x,»/,x / )|x-x'| 2 = J d V / (-A p )P T (0,r ] / ,p) lp=0 = B 2 T 1+ '' (54) 

If > defined in eq.(52) belongs to L\ oc (drj) then instead of the equalities 

in eqs.(53)-(54) we have bounds from above by Bit^ +v in cq.(53) and B 2 t 1+v 
in eq.(54). 

Proof:we prove eq.(54) (eq.(53) is simpler and proved in a similar way). Let 
us calculate 

(- A p ) / d V 'P r (0, rf, p)| p=0 = / d V >(2KT)-h cxp(-^(V) 2 ) 

E[t 1+ » $1 TrV (sr-ir?' + 7 («)) ds] = B 2 t 1+ » [ ' 

If W > then instead of the expectation value (55) we have 

fdr/(2^T)-iexp(-^(r/) 2 ) 



r 1+ - jl TrV (sr" V + 7 («)) dscxp (- / Q T (sr" 1 ?/ + V^ 7 (f )) da) 



< B 2 -, 
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where the inequality follows from W > 0. 

Corollary 3 means that if v < then the sample paths of diffusions generated 
by operators with singular coefficients have worse continuity properties than the 
Brownian paths (for Brownian paths see [8])). 

4 More general metrics 

We study the lower bound on G following from eq.(22) 
G L (0,0, P ) = J oo dr(2nr)-i 

exp (-t/q ds J dy{2irs(l - s))"i cxp(-2j^y)V (y/ryj) 
and the upper bound following from eq.(26) 

G^(0,0,p) 
=J™dr(2KT)-h 

Jo dsfdy(2irs(l - s))~i exp(- 2a( f_ 8) ) exp (-rV (y/ry)) 

for more general V 

A generalization of Theorem 1 reads 
Theorem 4 

Let us consider V = pVp which is not scale invariant but of the form 

V(v) - V{r,)f{ri) + 1( V ) (56) 

where V is a matrix scale invariant function (28) satisfying the conditions of 
Theorem 1 with — | < v < 0, / is a bounded function with a strictly positive 
lower bound, I is a nonnegative bounded matrix function. Assume in addition 
that 

J dyeM-\)f(y)V{y) >cl>0 (57) 

where c is a positive number. Under our assumptions (56)-(57) for any A > if 
|p| > A then there exist a positively definite bilinear form h% and constants K\ 
and K 2 such that 

K^phip)-" < G(0,0;p) < K 2 \p\- 2 " (58) 

If v > for V in eq.(56) then for |p| > A the inequalities (58) hold true with 
u = \. 

Proof: our assumptions (56) on V mean that it satisfies the inequalities 

T v pVt {y)p + php < V(V^y) < T v pV 2 {y)p + php (59) 
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with certain matrix functions V\ and V2 independent of r and bilinear forms l\ 
and h (independent of y) . It follows that the integral of P T satisfies the bounds 

J drP^ 2 cxp(-TpZ 2 p) < J drP T < J drP^ 1 exp(-rpZip) (60) 

where in the lower bound P L2 the potential V2 from the rhs of eq.(59) is applied 
and in P Ul the one from the lhs of cq.(59).The integral (57) defines h of cq.(32) 
(and the /i 2 from the upper bound (59)). Let us change the integration variable 

r = r{ph 2 p)^~ on the lhs of eq. (60) and r = p|p| _4 "|Vi(y)lo 2 " on the rhs. 
Then, the lower and upper bounds read (from eqs.(31),(34) and (38)) 

(ph 2 p)-^ dr(2irry? exp(-r 1+I/ - r{ph 2 p)- 2 "pl 2 p) 
<G(0 ) 0,p)<|p|- 2 «/ oo dp(27rp)-' 

Jo ds I du(2n)-Hs(l s))-^TW|o w exp(-£) ^ 61 ) 

cxp (- P i+» - pipr 4 "ifi( Uv ^(i^))io 2 "p^ip) 

The condition (57) implies that the bilinear form /i 2 is strictly positive. Hence, 
there exists a constant K such that 



Kph 2 p > php 

Then, for — \ < v < and |p| > A there exists c\ such that 

r(p/i 2 p) -T ^pl 2 p < rci 

in the exponential on the lhs of eq.(61). The l\ term can be set zero for the 
upper bound. In such a case for each A > there exist constants c\ and c 2 such 
that if |p| > A then the inequalities (61) take the form 



(phipY u J™ dr(27rr)-2 exp(-r 1+,y - rci) 

<G(0,0,p)<|p|- 2 - J^dp(2np)-i (62) 
J 1 d s jdu(2n)-H s (l- S ))-^\V 1 ( U )\^CM^)^V(-P 1+,J ) 

The inequalities (62) coincide with (58) because under our assumptions the 
integrals in cq.(62) are finite. The last statement of Theorem 4 follows from the 
inequalities (60) because the behavior for large p follows from the behavior of 
P T for a small r. If v > then in eq.(60) t 1+v < At for any A and a sufficiently 
small t. Hence, we obtain the same behavior of G for large momenta as in the 
case V = 1. 

We would like to note that the restrictive form (56) of V is not necessary. As 
an example we could consider V which has singularities at several points, e.g. 

V( V )=p 2 (a \ V - m \^+K\r 1 \ 2 n (63) 
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with 1 1/ 1 < M (only negative indices are non-trivial). An application of the 
lower and upper bounds (31) and (34) to the potential (63) leads to the conclu- 
sion that after an integration upon r the inequalities (40) hold true for |p| > A. 
Hence, the leading singularity v determines the behavior at large momenta. 

5 Discussion and summary 

As we pointed out in the Introduction our results concerning the Green func- 
tions can find an application to quantum field theory in an expanding universe. 
The stronger damping in momenta (eq.(40)) in the inflationary models {a > 1) 
at rj = 7]' = indicates that it would be promising to start quantization at this 
time (rj = corresponds to t = oo in cosmological models with a > |). The 
exponential inflation can be obtained as a limit a — > oo which corresponds to 
v = -§. This limit is beyond our rigorous approach but it could be treated by 
means of more sophisticated methods of ref.[9]. By a formal scaling argument 
we obtain again the behavior (40) which in inflationary cosmological models is 
known as the Harrison- Zeldovich spectrum of scalar fluctuations [1] [15]. The 
Green functions can be applied in order to derive a solution of Einstein equa- 
tions linearized around the homogeneous background [12] [13]. In such a case in 
addition to the scalar Green function the tensor Green function must be studied 
as well. Further consequences of our estimates concerning the spectrum of G 
for the complete theory still need to be explored. For this purpose a detailed 
dependence of the Green function on rj and rf would be useful. It is much harder 
to derive such estimates than the ones for the time zero case. In the Appendix 
we investigate the upper bound G u for general rj. In particular, calculations 
performed there suggest that it is only the behavior of V(y) for small y which is 
relevant for Theorem 4 and that the upper bound is valid for all v > — 1 + ^ . 
For the lower bound G L we can also obtain an integral representation. However, 
it is quite complicated. 

Another motivation for a study of the (d + l)-dimensional Green functions 
comes from the problem of a dimensional reduction of quantum fields defined on 
a brane [16]. In such a case we restrict ourselves to a (i-dimensional submanifold 
imposing the condition rj = rj' = 0. We have proved here that if the metric has 
a power-law behavior then the Green functions of the restricted quantum field 
theory are decaying faster in the momentum space than the standard |p| _1 . In 
particular, for v = — \ we obtain the propagator |p| -2 in d-dimensions which is 
the same as the one of the Euclidean massless free field. 

6 Appendix 

We calculate the upper bound for the Green's function G following form eqs.(26) 
and (13) in more detail. Set u — y/ry and let us perform the integration upon 



13 



t in eq.(26) with W = 0. Then 

G U (T),V',P) =2(2^-7/0 da Jdu(2ns(l-s))-*K ^2Mp7p (^ (64) 
where 

M = {ji — rj) 2 + (s (1 - a))" 1 (u - ?? - s (V - ??)) 2 (65) 

and K p is the modified Bcsscl function of order p [17]. 
The integral is simpler if 77 = 77' 

G^rj, r/, p) - 7T- 1 jl ds J duK (n^2 P V P (u^j (66) 

where 

N =(s(l- s))~*r}-u (67) 

If 77' = 77 = then the integral (64) further simplifies. If additionally V is scale 
invariant then we can calculate it exactly as in sec. 3. 

For arbitrary 77 and 77' the behavior of G u is much more complicated because 
the decay of G u substantially depends on 77. In the simplest case when V = \ 

G(r / ,77',p) = |pr 1 cxp(-|p||77'-77|) (68) 

Let us consider an arbitrary V and assume that it behaves as 

V{\y) = \ 2l 'B x {y) (69) 

when A — > with a certain B\ which as a function of A is bounded from above 
and from below,i.e.,C2(y) > B\{y) > C\(y) > for a small A. Let us change 
variables in eq.(64) 

U=\p\-Thy (70) 

Assume that |p| — > 00 , 77 — > and 77' — > in such a way that 9 = IpI 1 ^ 77 ?? and 
8' = |p|~77' remain finite. In such a case from eq.(64) we can conclude that 

G(\p\-rhe,\p\-rhe',p)^\p\-^ (71) 

for large p in agreement with eq.(40) and eq.(68) [y = for a constant V). 

Eqs.(64)-(67) give an integral representation of the upper bound which is 
expected to approximate the exact Green function G for large p. We suppose 
that the Fourier transform G u of G u is a reliable approximation to G at short 
distances. After the Fourier transform of eq.(26) with W = we can calculate 
the t integral in eq.(13) exactly. We obtain 

G^(77,x,77',x') = (27T)-** 2 fid8(s(l-s))-*fdy (detV (»))"* 

((x - x') V- 1 (y) (x - x') + (77 - T7') 2 + (s (1 - s))- 1 (y- v -a{rf- t?)) 2 ) ~ % 

(72) 
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If V % i(y) — 5^v(y) then the formula (72) can be expressed in a simpler form 

G u ( V ,x, 7?',x') = (27T)- 3 * 2 £ ( S (1 - / dy 

(jx - x'| 2 + v(y) (77 - i f + (s (1 - s))" 1 %) (y - 77 - s (77' - tj)) 2 ) 

(73) 

If u = 1 then eq.(73) gives 

(2 7 r)-^(|x-x'| 2 + (7 ? -77') 2 ) _V (74) 

as it should. 

The integrals (72)- (73) suggest some generalizations of the theorems proved 
in the main part. First, assume that V(y) ~ \y\ 2p for a large \y\ then the 
integrals (72)-(73) are finite (for |x — x'| 7^ 0) if p > — 1 + 3- Next, it can be 
shown from eq.(73) that if v{y) ~ \y\ 2u for y — > and 7/ = 77' = then 

G c/ (0,x, 0,x') ~ |x-x'r d+ TTT7 

for x — x' 0. The derivation of the result (44) based on eq.(73) suggests that 
for Theorem 4 only the behavior of V(y) for a small y is relevant (assuming the 
integral (73) is finite ). 

For general v and arbitrary x,x',7/ and 77' it is harder to obtain usable es- 
timates. Let us mention some special cases. It follows directly from eq.(41) 
that 

G(7?,x,0,x) ~ |, 7 |-<*(i+»0+i 

and 

G(0,x,t/,x)~ |77'r d(1+l/)+1 
whereas from eq.(73) we obtain that if v{rj) 7^ then 

G c/ (77,x,77,x')^7 J (77)-i|x-xr d+1 (75) 

when x — ► x'. 

If V(y) > c\y\ 2p with c > and p > — 1 + g for large ?/ then changing the 

integration variable y = |x — x'|~z we can show that for any 77 and 77' there 
exists A such that if |x — x'| > A then 

G u (r), x, 77', x') < if|x - x'r d+ ^ (76) 

When p > then eq.(76) gives a non-trivial estimate saying that the Green's 
function has a stronger decay for large distances than the one for operators with 
constant coefficients. However, such a decay at large distances will be changed 
by most perturbations W whereas the behavior for short distances is remarkably 
stable with respect to perturbations. 
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